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Proton-neutron, neutron- neutron and proton-proton momentum correlation lunctions (C pn , C„„, 
C pp ) are systematically investigated for 15 C and other C isotopes induced collisions at different 
entrance channel conditions within the framework of the isospin-dependent quantum molecular 
dynamics (IDQMD) model complemented by the CRAB (correlation after burner) computation 
code. 15 C is a prime exotic nucleus candidate due to the weakly bound valence neutron coupling 
with closed-neutron shell nucleus 14 C. In order to study density dependence of correlation function 
by removing the isospin effect, the initialized 15 C projectiles are sampled from two kinds of density 
distribution from RMF model, in which the valence neutron of 15 C is populated on both ld5/2 
and 2sl/2 states, respectively. The results show that the density distributions of valence neutron 
significantly influence nucleon-nucleon momentum correlation function at large impact parameter 
and high incident energy. The extended density distribution of valence neutron largely weakens 
the strength of correlation function. The size of emission source is extracted by fitting correlation 
function using Gaussian source method. The emission source size as well as the size of final state 
phase space is larger for projectiles sampling from more extended density distribution of valence 
neutron corresponding 2sl/2 state in RMF model. Therefore momentum correlation function can 
be considered as a potential valuable tool to diagnose the exotic nuclear structure such as skin and 
halo. 

PACS numbers: 25.60.-t, 25.70.Pq, 21.10.Gv 



I. INTRODUCTION 

Intensity interferometry method, developed by Han- 
bury Brown and Twiss (HBT) in 1950s [l|, was origi- 
nally used to measure astronomical objects such as an- 
gular diameter of stars. The method was later intro- 
duced into subatomic physics by Goldhaber et al., who 
extracted the spatial extent of an annihilation fireball in 
proton-antiproton reactions by two-pion correlations Q ■ 
Then the method was widely applied in exploring the 
nuclear reactions from low-energy to relativistic-energy 
[3l-fIol|. Recently, it has been extended to other fields, 
for instance, the analogous correlations describing the 
fermionic statics of elections |ll|, [l2[ . 

In heavy ion collisions (HICs) at intermediate energy, 
the HBT method is widely used to extract the space-time 
properties such as the source size and emission time of 
fragments by two-particle correlation functions [H, [TiJ ■ 
The dependences of momentum correlation function on 
the impact parameter [l5|, [l||, the total momentum of 
nucleon pairs [l|| [l?} , the isospin of the emitting source 
[l8j . the nuclear symmetry energy [l9| and the nuclear 
equation of state (EOS) [16| are also explored by exper- 
iment and theory. 

Besides the applications of two-particle momentum 
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correlation functions to investigate heavy ion collision 
process, the HBT method is also extended to study ex- 
otic structure of nuclei far from the /3-stability line due 
to the rapid development of radioactive nuclear beams. 
The neutron-neutron correlation functions (C, m ) of Bor- 
romean halo nuclei such as 6 He, n Li and 14 Be are con- 
structed to extract the size of separation between the 
two halo neutrons [20h25I |. In addition, Wei and Ma et 
al. found that the strength of proton-neutron momen- 
tum correlation functions at small relative momentum 
has linear dependence on binding energy per nucleon or 
sing le neutron separation energy for light isotope chains 
[lrj. I2IH [27l j . The analogous suppressed proton-proton 
correlation function is also suggested as another poten- 
tial tool to diagnose proton halo nuclei [28j beside the 
conventional methods such as total reaction cross sec- 
tion and momentum distribution width measurements 
[29T [30l | . Ma et al. have carried our experiment on 
RIPS in RIKEN to measure the proton-proton momen- 
tum correlation function for revealing the exotic struc- 
ture of proton-rich nucleus 23 Al |3ll. l32j. Therefore, it is 
very interesting to investigate how the exotic structure 
affects the nucleon-nucleon momentum correlation func- 
tion, which can serve as a potential observable to extract 
information about anomalous structure in nucleus. 

In this paper, we calculated 15 C + 12 C collisions by 
IDQMD model. 15 C is a one-neutron halo candidate be- 
cause of it's small neutron separation energy: S n — 1.218 
MeV [Hj], closed-neutron shell core: 14 C, narrower mo- 
mentum distribution of 14 C fragment from the breakup of 



15 C [34], |35[ and larger s-wave spectroscopic factor of 15 C 
ground state (3(| by 14 C(d,p) 15 C reaction measurement. 
However, a consistent picture have not been obtained 
in reaction cross section measurement. The interaction 
cross section (07) do not have peculiarity compared with 
neighbor isotopes at incident energy 740 MeV/nucleon 
[37J. However, reaction cross section (ctr) shows more or 
less enhancement at intermediate ener gies and there is 
also a large difference factor (d) for <jr |35l . l38j , which is 
defined as (39J: d — CT ^( c ^p)g^( G ) ; wnere an(exp) rep- 
resents intermediate energy experimental or and or{G) 
is the or calculated by the Glauber Model at the same 
bombarding energy with HO-type density distribution, 
which is obtained by fitting experimental or at high en- 
ergy. The s-wave spectroscopic factor extracted from 07 
is also different from the value obtained from transfer re- 
action. Fang et at found that s-wave component is domi- 
nant in the ground state of 15 C by simultaneous measure- 
ment of or and momentum distribution 35]. However, 
the fitted results by the Glauber model are deviated from 
the experimental data at low energies for or. Therefore, 
new probes are needed to estimate the density distri- 
bution of valence neutron in 15 C and other analogous 
neutron-rich exotic nucleus candidates. 

The initialized 15 C projectiles are sampled from den- 
sities calculated by RMF model to study how density 
distributions of outer neutron affect nucleon-nucleon mo- 
mentum correlation function. In RMF model, the last 
neutron of 15 C is populated on both ldh/2 and 2sl/2 
states, respectively. Because the isospin degree of free- 
dom is removed, the relationship between the momentum 
correlation function and the structure of exotic nuclei 
such as skin can be more directly explored by comparing 
different collisions induced by different configured 15 C 
projectiles. 

The rest part of this paper is organized as follows: in 
Sec. UH we briefly describe the HBT technique and the 
IDQMD model; the initialization of 15 C projectiles and 
nucleon-nucleon momentum correlation function of dif- 
ferent configured C isotopes induced collisions are dis- 
cussed in Sec. IIIII the summary is presented in Sec. IIVI 



tion g(p, x) in standard Koonin-Pratt equation 



II. HBT TECHNIQUE AND IDQMD MODEL 

A. HBT technique 

It is known that the final-state interaction (FSI) and 
quantum-statistical symmetry (QSS) affect the wave 
function of relative motion of light identical particles 
when they are emitted in close region in phase-space and 
time, which is the principle of intensity interferometry, 
i.e., the HBT method. The correlation function of two- 
particle can be obtained by convolution of emission func- 



C(P,q) 



Jd 4 Xl d 4 x 2 g(P/2,n) g(P/2,x 2 ) |0(q,r)| ; 
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where P(= pi + P2) and q[=( Pi - P2V2] are the to- 
tal and relative momentum of particle pair, respectively, 
g(p, x) is the probability of emitting a particle with mo- 
mentum p at space-time point x(r, t) and </>(q, r) is two- 
particle relative wave function with relative distance r = 
(r 2 - n) - |(vi + v 2 )(t a - h). 

In specific application of the Koonin-Pratt formula, the 
reliable single-particle phase space distribution at freeze- 
out stage is needed [15|. In this paper, the IDQMD model 
is used as event generator. It is a widely used transport 
model in intermediate energy HICs and has also been suc- 
cessfully applied to HBT studies for neutron-rich nuclei 
induced reactions by Wei and Ma et al. The 
phase space of emitted particles is used as the input of 
Pratt's CRAB (correlation after burner) code [io|, which 
takes the FSI and QSS effects into account for remedying 
the disadvantage of semi-classical transport model. 



B. IDQMD model 

The quantum molecular dynamics (QMD) approach is 
a many-body theory that can describe HICs from inter- 
mediate to relativistic energies (4l| - |43l |. The main ad- 
vantage of the QMD model is that it can explicitly treat 
many-body state of the collision system, so it contains 
correlation effects to all orders. Therefore, the QMD 
model provides valuable information about both the col- 
lision dynamics and the fragmentation process. The 
model also has excellent extensibility due to it's micro- 
scopic treatment of collision process. It mainly consists 
of several parts: initialization of the projectile and the 
target nucleons, propagation of nucleons in the effective 
potential, nucleon-nucleon (NN) collisions in a nuclear 
medium, the Pauli blocking and the numerical test. 

The IDQMD model is based on QMD model and af- 
filiates the isospin factors in mean field, two-body NN 
collisions and Pauli blocking [H, [H, 0-113 • In addi- 
tion, the phase space sampling of neutrons and protons 
in the initialization should be separately treated because 
of the large difference between neutron and proton den- 
sity distributions for nuclei far from the /3-stability line. 
In order to properly incorporate nuclear structure effects 
into microscopic simulations, stable initialized 15 C with 
and without neutron-halo structure have been sampled. 

In IDQMD model, each nucleon is represented by a 
Gaussian wave packet with width y/~L (here L = 2.16 
/to 2 ) centered around the mean position fl(t) and the 
mean momentum pl(t): 
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Then all nucleons interact via mean field and NN colli- 
sions. The nuclear mean field is parameterized by 

U( P , r z ) = a(^) + {3{^y< + hi r z )V c 
Po Po 2 



Po 

with po = 0.16 fm~ 3 (the normal nuclear matter den- 
sity), p, p n and p p are the total, neutron and proton 
densities, respectively. r z is zth component of the isospin 
degree of freedom, which equals 1 or -1 for neutrons or 
protons, respectively. The coefficients a, (3 and 7 are pa- 
rameters of nuclear equation of state (EOS). C sym is the 
symmetry energy strength due to the difference between 
neutron and proton, taking the value of 32 MeV. In this 
work, a = - 356 MeV, /3 = 303 MeV and 7 = 7/6 are 
taken, which corresponds to the so-called soft EOS with 
[42| . V c is the Coulomb potential and JJ Yuk is Yukawa 
(surface) potential, which has the following form: 

U Yuk = ^ V — exp(Lm 2 ) 
2m ^r? nj 

x [exp(— mrij)eiic(VLrn — rij/v&L) 

— exp(mrjj)erfc(vZm + r^/vlL)], (4) 

with V y = -0.0074 GeV, m = 1.25 /m" 1 and r v} = \n-fj\ 
is the relative distance between two nucleons. Experi- 
mental parameterized NN cross section, which is energy 
and isospin dependent, is used. 

The Pauli blocking effect in IDQMD model is also 
isospin dependent. The blocking of neutron and proton 
is separately treated as follows: each nucleon occupies 
a six-dimensional sphere with a volume of H 3 /2 in the 
phase space (considering the spin degree of freedom) and 
we calculate the phase space volume (V) of the scattered 
nucleons being occupied by the rest nucleons with the 
same isospin as that of the scattered ones. We then com- 
pare 2V /h 3 with a random number and decide whether 
the collision is blocked or not. 

The time evolution of the colliding system is given by 
the generalized variational principle. Nuclear clusters are 
constructed by a coalescence model, in which particles 
with relative momentum smaller than Po and relative 
distance smaller than Po are considered to belong to one 
cluster. The parameter set: Po = 300 MeV/c and Po = 
3.5 fm is taken here. 

III. THE INITIALIZATION OF 15 C 
PROJECTILES AND SYSTEMATICAL HBT 
RESULTS 

In IDQMD, the initialization of projectiles and tar- 
gets distinct proton from neutron. We sample the nu- 
cleon's coordinates from density distributions of proton 




R(fm) 

FIG. 1: (Color online) The proton, neutron, valence neutron and 
total density distributions of 15 C for the valence neutron on le£5/2 
and 2sl/2 states, respectively, by RMF calculations. 



TABLE I: RMS radii of initialized 15 C by RMF theory and 
experimental RMS charge radius of 14 C. The RMS radii of 
1 C projectiles are initialized to match with the values in this 
tabic. 
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and neutron, which are calculated by relativistic mean 
field (RMF) method. Fig. [T] shows the proton, neutron, 
valence neutron and total density distributions of two dif- 
ferent configured 15 C projectiles, whose valence neutron 
is assigned on both ld5/2 and 2sl/2 states, respectively. 
It can be seen that the valence neutron on 2sl/2 state 
has longer tail than that on ld5/2 state while proton 
almost has the same density distributions for two case. 
The effect of more extended valence neutron density dis- 
tribution can be reflected by nucleon-nucleon momentum 
correlation function as showed below. 

In our calculations, initialization of 15 C is carefully 
controlled. The stability of the sampled 15 C projectiles 
is strictly checked by time evolution in mean field till 200 
fm /c at zero temperature according to the average bind- 
ing energies, root-mean-square (RMS) radii and density 
distributions of neutron and proton. Eligible initializa- 
tion samples should meet the following requirements till 
200 fm/c: a) average binding energy needs to match 
with experimental data; b) RMS radius also needs to be 
in according with RMF result; c) the difference of neu- 
tron tails between two kinds of 15 C projectiles should 
keep all the way. To better imitate input density from 
RMF calculation and reflect the structure effect, thou- 
sands of eligible initialized samples are accumulated to 
simulating collisions. The initialized samples of other C 



4 




1 1 i 1 1 1 i 1 1 1 i 1 1 1 i 

(a) 

Incident energy: 800 Me V/u 
Impact parameter:7~8fm 

Y/Y rm >0 



C(narow outer neutron density) 
15 C(extended outer neutron density) 



40 60 80 
q(MeV/c) 



100 




40 60 

q(MeV/c) 



40 60 

q(MeV/c) 



100 



FIG. 2: (Color online) Proton-proton (C pp , panel (a)), neutron-neutron (C n n, panel (b)) and proton-neutron (C pn , panel (c)) momentum 
correlation functions as a function of relative momentum, which are calculated at incident energy of 800 MeV/nucleon, impact parameter 
6=7 — 8/m and selected nucleons with rapidity > 0. The circles and squares represent two different kinds of initialized 15 C projectiles 
induced collisions, where 15 C projectiles are sampled from two density outputs of RMF model with valence neutron on ld5/2 and 2sl/2 
states, respectively. The dashed line is used to guide the eyes. 
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FIG. 3: (Color online) (a): C pn strength at 5 MeV/c as a function 
of impact parameter, which is calculated at incident energy of 800 
MeV/nucleon and selected nucleons with rapidity > 0; (b): C' pn 
strength at 5 MeV/c as a function of incident energy, which is 
calculated at impact parameter 6 = 7 — 8/m and selected nucleons 
with rapidity > 0. The symbols in both panel (a) and (b) have the 
same convention as in Fig. [2] 



isotopes are prepared in a similar way above. 

Even with very elaborate initialization, the tail of den- 
sity in QMD model can not be reproduced very well com- 
pared with RMF model since the wave function in QMD 
model has Gaussian form. However, the obvious differ- 
ences of neutron tails between two kinds of 15 C projec- 
tiles can keep stable enough and thus it can play the 
role of skin structure in collisions. In the following, we 
will present how the skin structure expressed by different 
neutron density distributions in QMD affects nucleon- 
nucleon momentum correlations. 

The proton and neutron phase space of 15 C + 12 C col- 
lisions at freeze-out time generated by IDQMD model, 
is used as input to CRAB code. The obtained proton- 
proton, neutron-neutron and proton-neutron momen- 
tum correlation functions are shown in Fig. [2] respec- 
tively, where q denotes relative momentum of nucleon 
pair. As expected, C nn and C pn peak at small q while 
proton-proton is anti-correlation at small q owning to the 



Coulomb potential and antisymmetrization. The peak of 
C pp at 20 MeV/c is due to s-partial wave of the proton- 
proton scattering, which strongly depends on the size of 
emitting source. Our simulations pretty well reproduce 
the shape and height of C pp , C, 



and C pn vs. q com- 
pared with experimental cases |I4(. It can be seen that 
C nn and C pn are both largely reduced at small q because 
of the more extended neutron density distribution. In 
Fig. [21 the size of emission source is mostly decided by 
projectile-like remnants since we adopt impact parameter 
6 = 7 — 8/m and rapidity > 0. The RMS radii of proton, 
neutron, valence neutron, 14 C core, 15 C by RMF calcu- 
lation and experimental RMS charge radius of 14 C are 
shown in table [I] The RMS radii of initialized 15 C pro- 
jectiles are required to meet the RMF calculated values. 
Though there is no available experimental RMS radius 
of 15 C, RMS radius of 14 C core reproduces experimental 
RMS charge radius of 14 C quite well. The RMS radii 
of 15 C do not have large difference for two kinds of 15 C 
projectiles. However, the corresponding RMS radii of va- 
lence neutron are 3.85 fm and 5.01 fm. Therefore, the 
strength of momentum correlation functions can indeed 
sensitively reflect the fine difference of neutron density 
distribution. Although the density distributions of pro- 
ton for two configured 15 C are almost the same (as seen 
in Fig. [T]), Cpp still has large differences, which demon- 
strate that the proton density distribution has changed 
when proton emission. The density distributions of neu- 
trons and protons are coupling with each other during 
collisions. 

The impact parameter dependence of C pn at 800 
MeV/nucleon incident energy for the two different kinds 
of initialized 15 C projectiles is shown in panel (a) of Fig. 

The strength of correlation becomes decreasing with b 
up to 6.5 fm. In central collisions, emitted nucleons have 
stronger correlations between themselves because they 
come from one compact, hot and dense region, which is 
consistent with previous Boltzmann-Uehling-Uhlenbeck 
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FIG. 4: (Color online) The same as Fig. [3] but for C n 



FIG. 5: (Color online) The same as Fig. [3] but for C p 



(BUU) PJ, Eg and IDQMD [JJ] results. In peripheral 
collisions with b > 6.5 fm the rapid enhancement of C pn 
reveals different collisions dynamics from central colli- 
sions. This is due to the fact that there is no bulk overlap 
region at large impact parameter and only the outer nu- 
cleons of 15 C are scraped. They keep more initial struc- 
ture information of projectile. Thus, there exists strong 
correlation among these nucleons. The difference of 15 C 
in neutron density distribution such as skin structure can 
be well revealed in peripheral collisions while the small 
difference of density is wiped away in violent central and 
semi-peripheral collisions. 

How incident energy affects C pn is represented in panel 
(b) of Fig. [3] The correlations for both cases increase 
from very low incident energy and then almost reach sat- 
uration above 300 MeV/nucleon. This is understandable 
due to the fact that increasing incident energies lead to 
more rapid collision process and there is smaller space 
and time interval among emitting nucleons [THIHI, which 
leads to stronger correlation. The apparent differences 
in C pn exist for the two kinds of 15 C induced collisions 
above saturation energy. The saturation of C pn at high 
incident energies and the large differences provide us the 
proper entrance channel conditions to explore the rela- 
tion between initial structure and final effect. 

The impact parameter and incident energy depen- 
dences of C' nn are shown in panel (a) and panel (b) of 
Fig. HI respectively. The tendency of C nn is similar to 
Cp n for both impact parameter and energy dependences. 
The analogous results are also obtained for C pp shown in 
Fig. El 

The strength of C pn , C nn and C pp is mainly decided by 
the distance in phase space at freeze-out time. Under spe- 
cific entrance channel conditions (high energy and large 
impact parameter), the size of phase space in final state is 
consistent with initial size of 15 C projectile, which largely 
depends on the neutron skin structure. This can be seen 
from Fig. E] showing impact parameter dependences of 
RMS radii in coordinate and momentum space. The size 
of phase space decreases with the increasing of impact pa- 
rameter. The more extended neutron skin leads to larger 
RMS radii of coordinate and momentum in final state 
for peripheral collisions at 800 MeV/nucleon, which re- 
sults in smaller correlation in the HBT method. However, 
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FIG. 6: (Color online) Root-mean-square of coordinate (panel (a) 
and (c)) and momentum (panel (b) and (d)) space distributions in 
the center of mass reference frame as a function of impact param- 
eter at incident energy of 800 MeV/nucleon and rapidity > for 
neutron (panel (a) and (b)) and proton (panel (c) and (d)), respec- 
tively. The solid (empty) circles, squares and triangles represent 
transverse, parallel and total RMS radii, respectively, for 15 C pro- 
jectiles induced collisions, where 15 C projectiles are sampled from 
two density outputs of RMF model with valence neutron on ld5/2 
and 2s 1/2 states, respectively. 



RMS radii do not have difference between the two kinds 
of 15 C induced reactions in central and semi-peripheral 
collisions because the only difference in outer neutron is 
covered in violent collisions. Therefore, the physical pic- 
ture of phase space size is consistent with impact param- 
eter dependences of C prl , C nn and C pp . The excitation 
function of RMS radii of phase space in peripheral col- 
lisions shown in Fig. [7] gradually become saturated at 
high incident energy, which can explain the saturation of 
correlations in panel (b) of Fig. El 2] and El 

The space-time extent of emission source for different 
particles can be extracted from the shape and height of 
Cpp, C nn and C pn for 15 C projectiles with and without 
neutron skin. Angle-averaged Koonin-Pratt formula can 
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FIG. 7: (Color online) Root-mean-square of coordinate (panel (a) 
and (c)) and momentum (panel (b) and (d)) space distributions in 
the center of mass reference frame as a function of incident energy 
with impact parameter b = 7 — 8/m and rapidity > for neutron 
(panel (a) and (b)) and proton (panel (c) and (d)), respectively. 
The symbols have the same convention as in Fig. \E\ 
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FIG. 9: (Color online) Dependence of C pp (strength of C pp at 20 
MeV/c) on C isotopes ( 12 - 13 > 14 ' 15 > 16 C) induced collisions. Differ- 
ent collision systems are compared with the same reduced impact 
parameter range: 0.875 - f.O and only protons with rapidity > 
are used to calculate momentum correlation function. Nucleons 
Q f 12,13,14,16(2; are sam pl e( j f r0 m uniform density. We construct 
four kinds of 15 C projectiles: two of them are sampled from two 
RMF densities shown in Fig. [l] respectively, and the other two are 
sampled from uniform distribution with two different sizes, corre- 
sponding to the sizes of ldb/2 and 2sl/2 states in RMF model, 
respectively. 
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FIG. 8: (Color online) Gaussian source fitting for C pp (panel (a)), 
Cnn (panel (b)), C pn (panel (d)) and variance between IDQMD + 
CRAB correlation and Gaussian source correlation as a function 
of Gaussian source size (panel (d)). The entrance channel condi- 
tions of 800 MeV/nucleon bombarding energy and impact param- 
eter with b = 7 — 8/m are used and the nucleons with rapidity > 
in final state are chosen. 



be written as [H [5(| : 

R(q)=4n [ r 2 drK(q,r)S{r), 



(5) 



where S(r) is the isotropic source function standing for 
the probability distribution of emitting a pair of nucleons 
at relative distance r. K(q,r) is angle-averaged kernel 
function obtained from radial part of two-nucleon relative 
wave function. If we assume S(r) has simple Gaussian 
form: 



S(r) 



(27r) 3 / 2 a 3 



exp(- 



2cr 2 



(6) 



where a describes the spatial distribution of nucleons 
emitting location and the RMS radius of Gaussian source 
equals \/3a. Gaussian source method is used to fit the 
correlation functions (C pp , C nn and C pn ) for two kinds of 
15 C induced collisions. The panel (d) of Fig. [5]represents 
the variance between IDQMD and Gaussian source cor- 
relations as a function of RMS radius of Gaussian source. 
Then the best fitting is judged by the variance and the 
RMS radii obtained are shown in panel (a), (b) and (c), 
respectively. We get the same Gaussian source size from 
C nn and Cpn- The emission source size of projectiles 
with skin structure is about 20% larger than that of pro- 
jectiles without skin structure for C nn and C pn . The 
source size of C pp is smaller than that of C nn and C pn 
for both kinds of 15 C projectiles. Comparing with ini- 
tialized neutron distribution shown in Fig. [I] and table HI 
we can conclude that the initial structure information of 
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neutron skin and size of projectile is kept till final state. 
The size of proton emission source is different for the two 
kinds of 15 C projectiles, though initial density distribu- 
tion of proton is the same. The size is about 15% larger 
for projectiles with neutron skin, which proves that pro- 
ton density distribution is disturbed through interaction 
with neutron even in very peripheral collisions. 

The size of emission source extracted above is under 
hypothesis that nucleon is expressed by point particle. 
However, in QMD model, nucleon wave function has 
Gaussian form in coordinate and momentum space. The 
true distribution of source is the convolution of the distri- 
bution of wave packet center with Gaussian density dis- 
tribution of single nucleon. In Ref. [HI HU , the relation 
between RMS radius of the freeze-out points < r{t) 2 y 1 / 2 
and the variance of the source and wave function is ob- 
tained by assuming that a chaotic Gaussian source has 
formed in collision without correlation between coordi- 
nate and momentum space: 

< r(i) > 1 / 2 = y/3(A + L), (7) 

where A is variance of chaotic emission source with Gaus- 
sian form and L is square of width of Gaussian wave 
packet in QMD model with value of 2.16 fm 2 in our cal- 
culations. Now with finite width modification of Gaus- 
sian wave packet, then we get smaller RMS radius of 
emission source: V3A — ^/3(er 2 — L). For example, 
RMS radii of C pp are 5.12 fm and 6.07 fm for 15 C projec- 
tiles with and without neutron skin, respectively. RMS 
radii of both C nn and C nn have the same value: 5.88 fm 
and 7.37 fm, respectively. 

We also study the systematical dependence of mo- 
mentum correlation on other C isotopes induced colli- 
sions. Fig. [9] shows the strength of C pp at 20 MeV/c 
£ or 12, 13,14,15,16^ _|_ 12q S y S t ems . xhe collisions are com- 
pared at the same reduced parameter range: 0.875 - 1.0. 
Four kinds of 15 C projectile are constructed. The other 
two kinds of 15 C projectiles are sampled from uniform 
density distributions beside the sampling from RMF dis- 
cussed above. Nucleons are also uniform distribution in 
12 C, 13 C, 14 C and 16 C projectiles. The total systematical 
tendency of mass dependence of C pp is decreasing except 
for projectiles with neutron skin structure, which can be 
interpreted as following: though they are simulated at 
the same reduced impact parameter, 12 C + 12 C collisions 
have smaller emission source as compared with 16 C + 12 C 
collisions. Correlation for projectiles with neutron skin is 
stronger than other cases. This is because that the outer 
density distribution is more extended due to neutron skin 



and then it results in smaller overlap zone between pro- 
jectile and target for the same reduced impact parame- 
ter at peripheral collisions. Therefore, proton-proton has 
stronger correlation due to they come from one compact 
projectile-like remnants, which is obvious different from 
the case without skin structure. 

IV. SUMMARY 

The very small nucleon separation energy of last nu- 
cleon and suddenly increased radius compared with its 
neighboring isotopes are two main features of exotic halo 
nucleus. HBT method has been used to study the rela- 
tionship between momentum correlation and single nu- 
cleon separation energy several years before. In this pa- 
per, we systematically investigate how the density distri- 
bution of valence neutron in one-neutron halo candidate 
15 C affects the strength of nucleon- nucleon momentum 
correlation function. Specifically, two kinds of 15 C pro- 
jectiles are sampled from two different density outputs of 
RMF model, respectively. The difference of outer neu- 
tron density distribution between two kinds of initialized 
15 C samples can be viewed as neutron skin structure in 
IDQMD simulation. The more extended density distri- 
bution due to outer neutron in projectile results in the 
larger and hot emission source, which leads to a larger 
size of phase space in final state and then corresponds 
to weaker correlation. The energy and impact parame- 
ter dependences show that nucleon-nucleon momentum 
correlation function is a very sensitive observable of den- 
sity distribution of valence neutron at high energy and 
peripheral collisions. Therefore, momentum correlation 
function at high bombarding energy and large impact pa- 
rameter can serve as a new potential probe to diagnose 
the exotic structure such as skin and halo besides the 
traditional measurements of total reaction cross section 
and momentum distribution. 
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